Inhibited spontaneous emission of quantum dots observed in a 3D photonic band gap 
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We present time-resolved emission experiments of semiconductor quantum dots in silicon 3D 
inverse-woodpile photonic band gap crystals. A systematic study is made of crystals with a range 
of pore radii to tune the band gap relative to the emission frequency. The decay rates averaged 
over all dipole orientations are inhibited by up to a factor of 12 x in the photonic band gap, and 
enhanced up to 2x outside the gap, in quantitative agreement with theory. We discuss the effects 
of spatial inhomogeneity, non-zero non-radiative decay, and transition dipole orientations on the 
observed inhibition in the band gap. 



In the field of cavity quantum electrodynamics (QED) 
it has been recognized that the nano-environment of a 
two-level quantum system may serve to tailor the fun- 
damental light-matter interactions P, @] ■ Of particular 
interest is the broadband and radical suppression of vac- 
uum fluctuations in a three-dimensional (3D) photonic 
band gap, i.e., a frequency range for which light is for- 
bidden for all wave vectors and all polarizations (3j . Such 
band gaps are expected in 3D photonic crystals, i.e., di- 
electric nanostructures with periodicities less than half 
the optical wavelength [3]. Anticipated cavity QED ef- 
fects of band gaps include complete inhibition of spon- 
taneous emission, photon-atom bound states, collective 
laser-like emission, and intriguing fractional decay 

To date, the cavity QED effects of 3D photonic band 
gaps on two-level light sources have only been studied 
in theory 0-0]. In these studies, one usually considers 
a single two-level source (or N identical ones) with ideal 
100% quantum efficiency that is excited once at t=0, and 
that is placed in a perfect and infinitely extended pho- 
tonic crystal. Clearly, this idealized situation differs from 
real situations. Real sources are inhomogeneously broad- 
ened, have non-zero non-radiative decay, are repeatedly 
excited in time, and are embedded in photonic crystals 
of finite extent. In this paper, we perform a first exper- 
imental study of cavity QED effects of real 3D photonic 
band gap crystals on embedded two-level light sources. 
To this end, we have developed Si photonic crystals with 
a diamond- like structure that have broad gaps [l3|. We 
perform time-correlated single photon counting experi- 
ments on embedded semiconductor quantum dots and 
observe broadband inhibition inside and enhancement 
outside the band gap. We notably find that the effects of 
spatial inhomogeneity, non-zero non-radiative decay, and 
transition dipole orientations are relevant in the band 
gap. Our results allow for the first time to test predic- 
tions from several different theoretical models. 

Silicon 3D inverse woodpile photonic band gap crys- 
tals were fabricated by a CMOS-compatible method de- 
scribed in Ref. [Hj . In brief, two orthogonal sets of pores 
are etched consecutively in a silicon wafer by reactive ion 
etching after careful alignment. Fig. [Tfa) shows a scan- 



ning electron micrograph of a typical crystal. The cubic 
crystals have lattice parameters a = 693 and c = 488 nm 
(a/c = and a range of pore radn (136 < r < 186 
nm) to tune the band gap relative to the emission spec- 
trum. The 3D crystal extends over L"^ = 12 x 12 x 12 ac^, 
which exceeds the Bragg length in every direction 
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The good optical quality and high photonic strength of 
our crystals have been confirmed by optical reflectivity 
where intense and broad peaks are seen. The stopbands 
overlap for all probed directions and polarizations, which 
is an experimental signature of a photonic band gap 12 j . 

Fig. [Hb) shows the calculated band structure [l^ for 
one of our inverse woodpile crystals. A broad band gap 
appears in the frequency range where modes are forbid- 
den for all wavevectors. Fig. [TJc) shows the correspond- 
ing density of states (DOS). At low frequency the DOS in- 
creases quadratically as is typical of an effective medium. 
Beyond 0.35 the photonic regime is entered as the DOS 
deviates from the parabola. The DOS completely van- 
ishes in the photonic band gap between 0.478 and 0.504. 

To study emission from two-level sources, we immersed 
the crystals in suspensions of PbS colloidal quantum dots 
in toluene (Evident). The concentration was as low as 
2 • 10~^ M to avoid dot-dot interactions such as energy 
transfer. The dots emit at photon energies between 0.8 
and 0.9 eV including the telecom range (Fig. 0]), and 
their transition dipole orientations sample all directions. 
We measured room temperature decay rates at three en- 
ergies in the quantum dot spectrum: 0.828, 0.850, and 
0.893 eV (A = 1500, 1460, 1390 nm respectively). Out- 
side the crystals the dots reveal single exponential decay 
with an energy-independent rate. By keeping the dots 
in suspension the photophysical properties remain stable 
for months, much longer than the short stability of dried 
dots in previous studies [l^ . 

The quantum dots were excited with short (11 ps) light 
pulses from a frequency-doubled Nd"^+ : YAG laser (Time 
Bandwidth Cougar) at A = 532 nm at a repetition rate 
of 409 kHz. The laser power was kept sufficiently low 
to avoid saturation and ensure that the dots remain in 
the linear regime. Light collected from the quantum dots 
was resolved by a monochromator (Acton 25001) set to 
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Figure 1: (color) (a) Scanning electron micrograph of a 3D 
inverse woodpile photonic crystal made from silicon. The 
crystal consists of two perpendicularly etched sets of care- 
fully aligned pores, and is delimited by the dashed lines, (b) 
Band structures of an infinite inverse woodpile photonic crys- 
tal with pore radius 170 nm, using e = 12.1 for silicon and 
e = 2.25 for toluene-filled pores. The band gap is indicated 
with the red bar. Inset: first Brillouin zone, (c) Density 
of states (DOS) per volume for the same crystal calculated 
with 2000 fc-points [l3|. The DOS vanishes in the band gap. 
Dashed curve: quadratic behavior in the low frequency limit. 



0.003 eV bandwidth, sufficiently narrow to resolve the 
photonic gaps. The excitation light was focused with 
an NA=0.12; 4x objective on the crystal, see Fig. [21[a). 
To minimize background from the suspension, light from 
the quantum dots was collected at 90° from the excita- 
tion beam by a detection objective with NA=0.7; lOOx. 
To precisely align the crystal the excitation focus was 
scanned along the edge of the structure. Fig. ^h) shows 
a scan of intensity versus y-position of the excitation fo- 
cus. We see a sharp transition from high (in suspension) 
to lower intensity in the crystal. In a scan performed on 
bulk silicon next to the crystal a sharp edge is also seen. 
Here the intensity is about half of the intensity collected 
from the photonic crystal. Hence about half of the inten- 
sity is emitted by quantum dots in the photonic crystal 
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Figure 2: (a) Drawing of the crystal position relative to exci- 
tation and and detection foci. The coordinate system of the 
position scans is shown, (b) Peak intensity of quantum dots 
at 0.842 eV (1475 nm) vs. y-position in an inverse woodpile 
crystal (red squares) and near a silicon wafer (black circles), 
curves are guides to the eye. The photonic crystal reveals 
twice as much signal as bulk silicon. The measured intensi- 
ties {PC, S) are input for the time-resolved emission model. 



while the other half stems from quantum dots in suspen- 
sion. Therefore, we conclude that a significant sponta- 
neous emission signal originates from quantum dots in 
our photonic band gap crystals. 

To precisely measure the arrival time of the emitted 
photons, time-correlated single photon counting was em- 
ployed with a cooled photomultiplier tube (Hamamatsu 
NIR) and a timing card (Picoquant PicoHarp 300). Since 
we perform experiments at relatively low photon energy 
in the near infrared to avoid Si absorption, the detec- 
tor darkcount is substantial (3 • 10*^ counts/s) and the 
quantum dot signal is relatively low (< 1% of the repeti- 
tion frequency, or < 4-10'^ counts/s). Therefore, we have 
collected for long times (hours) to obtain sufficient statis- 
tics, and carefully subtracted the background [15|. Fig. [3] 
shows time-resolved spontaneous emission for quantum 
dots in two different photonic crystals with pore radii 170 
and 136 nm. The emission inside the bandgap (r = 170 
nm) decays slower than outside the bandgap (r = 136 
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Figure 3: Time-resolved spontaneous emission measured at 
0.850 eV on two different photonic crystals with pore radii 
170 nm (closed squares) and 136 nm (open circles). Bi- 
exponential models (Eq. [TJ are indicated by the curves (solid 
for 170 nm, dash-dotted for 136 nm). Bottom panel: the 
residuals are random indicating high quality fits. 



nm), which confirms that the excited state lifetime of 
quantum dots is controlled by the photonic crystals. 

Interestingly, Fig. [3] reveals that the photon arrival 
times decay monotonously in time to a vanishing inten- 
sity, in agreement with Ref. Q; we do not detect sig- 
natures of fractional decay or oscillations predicted in 
Refs. jTi^tS]- A possible reason for this discrepancy may 
be that those predicted features are not robust to ensem- 
ble averaging. Other reasons could be the finite dynamic 
range (1.5 decades) in our experiments, the assumption 
of very large oscillator strengths in some models, or that 
theories predict excited-state population dynamics in- 
stead of photon arrival times while these two phenomena 
do not necessarily have the same dynamics, see Ref. flij . 

Since the emission originates both from quantum dots 
in the crystal and in suspension, we model the time- 
resolved emission with a double-exponential: 

fit) = I{SeM-lst) + PCexpi—fpct)). (1) 

Here jpc is the rate of quantum dots in the photonic 
crystal that we wish to obtain, and 7s is the rate of 
dots outside the crystal that is separately found to be 
7s = 1.9 /is^^. S and PC are determined by scans of 
the detection focus (see Fig. [2lb)). The high quality of 
our model is confirmed by residuals randomly distributed 
about (Fig. [3] bottom) , and by the goodness of fit xled 
of 1.01 and 0.93 near unity. From Fig. [3] we obtain de- 
cay rates 7pc = 3.1 ^s~^ and 0.88 fis~^ for two dif- 
ferent crystals. Therefore these data reveal respectively 
enhanced and inhibited emission compared to quantum 
dots in suspension. 

Fig. U (a) to (c) collects the decay rates measured on 
several crystals with a range of pore radii and hence gap 
frequencies. To investigate the inhibition in the photonic 
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Figure 4: (a) to (c): measured inhibition 'ys/lPC (open cir- 
cles) versus photon energy for crystals with indicated pore 
radii. Each vertical bar spans variations of 2 to 3 measure- 
ments taken on a same crystal on different days. The horizon- 
tal lines indicate reference level 1. Circles represent extremal 
data since a lower inhibition is likely biased by dots outside 
the crystal. The calculated inverse density of states (DOS) 
are shown as curves. The top panel shows the quantum dot 
emission spectrum. 



band gap range we plot the ratio "fs/lPC- For all crys- 
tals a good agreement between the measured inhibition 
and the inverse DOS is observed. The agreement is much 
better than previously found in inverse opals [13, likely 
since our quantum dots occupy a much larger fraction of 
the unit cell thus providin g a better average of the local 
density of states (LDOS) [l3|- For the 145 nm crystal 
(Fig. Ufa)) the photon energy is higher than the band 
gap frequency. An enhanced decay rate up to 2x is ob- 
served with increasing photon energy, in agreement with 
the DOS. The 170 nm crystal (Fig.gljb)) reveals up to a 
2-fold inhibition at frequencies just below the gap. Strik- 
ingly, a strong up to 12 x inhibition is found at photon 
energies deep in the band gap where a minimum decay 
rate of 0.16 ns~^ is observed. The corresponding life- 
times up to Ti = 6.25 ^s are extremely long for quantum 
dots, and confirm strongly shielded vacuum fluctuations, 
which is favorable for applications in quantum informa- 
tion processing [l^l ■ The occurrence of a lower inhibition 
in the band gap may be caused by slight misalignments 
that bias results to 7s; since the measurements take long 
it is conceivable that the alignment slowly drifts during 
data collection. The 176 nm crystal (Fig. |H[c)) reveals 
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Figure 5: (color online) Envelope of the LDOS normalized to 
vacuum versus depth into a finite band gap crystal with extent 
L = 12a. In the band gap the LDOS decreases exponentially 
away from the crystal surface to a minimum < 10~^, as es- 
timated from Ref. [l^l- Outside the gap the LDOS remains 
finite, as shown for a pseudogap as in Ref. [13 ]. 



a clear inhibition of the decay rate inside the band gap 
in good agreement with the calculated density of states. 
Up to 12x inhibition is observed, which agrees with the 
results on the 170 nm crystal. 

All theoretical studies predict strongly inhibited 
excited-state lifetimes in the range of the band gap 
[^. Indeed we observe drastically inhibited emission in 
our photonic band gap crystals. Interestingly, the sit- 
uation studied here leads to new physics since emission 
rates in real 3D photonic band gaps will reveal spatial 
inhomogeneity. Based on theory for 2D crystals 19] we 
expect the LDOS averaged over the unit cell to decrease 
exponentially with depth z in the crystal at frequencies 
in the band gap, see Fig. [SJ In contrast, outside the gap 
the LDOS varies only slightly with z from vacuum to the 
finite bulk value (c/. Fig. HJc)). In the gap, the LDOS 
is minimal at the crystal center with a minimum deter- 
mined by the crystal size L [2C|. Hence we estimate the 
maximum inhibition at the center of our crystals to be 
more than 100 x, see Fig. [5] 

Spatial inhomogeneity has intriguing consequences: 
experimentally we measure time-resolved emission that is 
proportional to the radiative rate ^rad times the LDOS- 
distribution [16[. Thus such an experiment is biased to 



quantum dots with a high emission rate and against dots 
with a strong inhibition. This qualitatively explains why 
we observe a finite inhibition. To take advantage of the 
strong inhibition deep inside 3D crystals calls for com- 
pletely different experimental approaches, viz., the prob- 
ing of excited-state populations to directly verify the the- 
oretical predictions [H-Q. We propose to use transient 
absorption as a probe of excited states that become more 
and more stable due to the shielding of the vacuum fluc- 
tuations in a 3D photonic band gap. 

In addition to spatial inhomogeneity, non-zero non- 



radiative decay occurs for any real light source. Hence, a 
time-resolved emission experiment yields a total rate that 
equals the sum of the radiative and non-radiative rates: 
Itot = Irad + l7irad- When ^rad IS strougly inhibited in a 
3D photonic band gap, as is the case here, the rate ^tot is 
minimal with a lower bound given by jnrad- While jnrad 
is not well known for our quantum dots, we can estimate 
an upper bound from the maximum observed inhibition 
to be 1/12 X 1.9 — 0.16/is^^, which is reasonable in view 
of data on other colloidal dots [2l| . 

Finally, it is well-known that the emission rate strongly 
depends on orientation 22] . Since our quantum dots are 
in suspension they sample all orientations of the transi- 
tion dipole moment. Thus the inhibition observed here 
is robust to orientation averaging. This is in contrast 
to recently observed inhibited emission in 2D photonic 
crystals [isj, and nanowires f23| that concern dots with 
a particular orientation of the transition dipole moment. 
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